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Abstract 

This  paper  develops  a  new  thermodynamic  model  that  predicts  maximum  fuel-cell  efficiency  and  fuel  utilization  as  a  function  of  fuel  composition 
and  operating  conditions.  Interestingly,  it  is  shown  that  maximum  possible  efficiency  is  independent  of  membrane-electrode  assembly  characteristics 
and  internal  polarization  losses.  However,  because  the  power  density  depends  on  internal  polarization  losses,  the  cell  size  needed  to  achieve  maximum 
efficiency  depends  greatly  on  cell  architecture.  A  very  detailed  fuel-cell  model  is  used  to  validate  the  thermodynamic  model  and  to  provide  further 
insight  on  practical  considerations  such  as  power  density.  The  model  is  illustrated  using  solid-oxide  fuel  cell  examples. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

It  is  well  known  that  fuel  cells  promise  the  potential  for  very 
high  efficiency  in  the  conversion  of  hydrogen  or  hydrocarbon  to 
electricity.  However,  precise  definitions  of  efficiency  can  vary 
widely  [1-3].  In  a  narrow  sense,  an  efficiency  can  be  defined  for 
the  fuel  cell  alone.  In  a  broad  sense,  efficiency  must  consider 
the  entire  system,  including  all  balance-of-plant  components. 
Because  full-system  designs  can  vary  greatly,  it  is  difficult,  if  not 
impossible,  to  generalize  the  analysis  of  efficiency.  The  objective 
of  the  present  paper  is  to  develop  purely  general  thermodynamic 
models  that  characterize  efficiency  of  the  fuel  cell  alone. 

An  interesting  and  remarkable  result  of  the  thermodynamic 
analysis  is  that  maximum  fuel-cell  efficiency  is  independent  of 
internal  polarization  losses.  However,  power  densities  depend 
strongly  on  internal  losses.  To  achieve  maximum  efficiency,  the 
cell  size  depends  on  power  density  and  hence  on  the  struc¬ 
ture  and  performance-characteristics  of  the  cell.  To  assist  in¬ 
terpretation  of  the  thermodynamic  models,  this  paper  also  uses 
a  very  detailed  model  of  solid-oxide  fuel  cells  (SOFC).  The  de¬ 
tailed  model  considers  specific  membrane-electrode  assembly 
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(MEA)  structure,  reactive  porous-media  transport,  electrochem¬ 
ical  charge  transfer,  and  elementary  heterogeneous  reforming 
chemistry  [4]. 

The  thermodynamic  analysis  developed  in  this  paper  is  both 
general  and  restrictive.  It  is  general  in  the  sense  that  it  accom¬ 
modates  different  fuels,  flow  rates,  temperatures,  and  pressures. 
Because  maximum  efficiency  and  fuel  utilization  are  thermody¬ 
namic  functions,  they  do  not  depend  on  particular  MEA  structure 
or  cell  design.  However,  the  analysis  is  also  restrictive.  Isother¬ 
mal  operation  is  assumed  and  all  gases  entering  the  fuel  cell  are 
assumed  to  be  at  the  cell  operating  temperature.  Thus,  the  cost  of 
heating  inlet  streams  (especially  the  air)  must  be  accommodated 
in  balance-of-plant  models  that  consider  heat-exchanger  perfor¬ 
mance.  The  models  presume  that  charge  transfer  proceeds  via 
Ht.  When  the  parent  fuel  is  a  hydrocarbon  internal  reforming  is 
assumed  to  proceed  to  equilibrium. 

Although  the  thermodynamic  model  has  restrictions  and  as¬ 
sumptions,  the  results  are  highly  informative  and  practically  use¬ 
ful.  Overall  thermal  management  of  an  entire  system  is  critically 
important.  However,  the  thermodynamic  model  developed  here 
can  be  integrated  into  a  larger  system  model.  The  assumption  of 
hydrogen  charge  transfer  is  a  very  reasonable  one,  and  does  not 
affect  the  overall  greatly. 

The  thermal  efficiency  of  any  energy-conversion  device  can 
be  defined  as  the  ratio  of  useful  work  done  Wout  and  the  potential 
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of  the  inlet  stream  to  do  work  Qm\  that  is  e  =  Woui/ (QAn.  For  a 
fuel-cell  system,  the  useful  work  is  electric  power  We,  which  is 
the  product  of  the  electric  current  density  i  and  operating  voltage 
fc’cdi,  integrated  over  the  active  ME  A  area,  We  =  f  iEceudA. 
The  potential  to  do  work  can  be  represented  in  terms  of  the  heat- 
release  associated  with  full  oxidation  of  the  inlet  fuel  stream, 
<2in  =  wf  in  A/i[jn.  Here  the  inlet  fuel  mass  flow  rate  is  rhfjn 
and  A/if  jn  is  the  specific  enthalpy  associated  with  completely 
oxidizing  the  fuel  stream.  With  these  definitions  in  mind,  the  net 
fuel-cell  efficiency  can  be  defined  as, 

We  f  i Ziceii  d A 

8=  —  =  ^ - — - .  (1) 

Q  in  ^f,in  A^f,in 

The  efficiency  of  a  fuel-cell  system  can  be  written  as  the  product 
of  three  contributing  efficiencies  [5  ] :  the  reversible  efficiency  £r, 
the  voltage  or  part-load  efficiency  gy,  and  the  fuel  utilization  £y, 

S  =  £r£V£U-  (2) 


The  ideal  efficiency,  or  the  reversible  efficiency,  £r  is  written  as 


AG  AS 

£r  =  -  =  1  -  T - , 

AH  AH 


(3) 


where  AG,  AH,  and  AS  are  the  changes  in  molar  free  energy, 
enthalpy,  and  entropy,  respectively,  associated  with  full  oxida¬ 
tion  of  the  fuel. 

Any  overpotential  losses  within  the  fuel  cell  reduce  the  cell 
potential  when  it  is  operated  under  load.  Thus,  the  net  efficiency 
depends  on  the  operating  cell  potential  Eceii.  A  part-load  effi¬ 
ciency  or  the  voltage  efficiency  is  defined  as 


£V  = 


£cell 

Erev 


(4) 


where  Eiev  is  the  reversible  (Nernst)  potential.  Thorstensen  [5] 
proposed  that  overall  cell  efficiency  should  be  the  product  of  the 
reversible  efficiency  and  the  part-load  efficiency  as. 


£  =  £r  £y- 


(5) 


Ferguson  and  Ugursal  improved  this  approach  by  incorporating 
activation-polarization  losses  using  a  two-linear-segment  polar¬ 
ization  curve  [6], 

Fuel  utilization  also  strongly  affects  the  efficiency.  Consider 
an  SOFC  system,  where  fuel  is  electrochemically  oxidized  along 
the  length  of  an  anode  channel.  As  the  fuel  is  consumed,  the 
anode  fuel  stream  is  also  diluted  by  reaction  products  (i.e.,  H2O 
and  CO2).  As  the  fuel  concentration  decreases  along  the  length 
of  the  anode  channel,  the  reversible  potential  £rev  decreases. 

As  long  as  the  reversible  potential  exceeds  the  cell  operating 
potential  (i.e.,  Eiev  >  /iceii  )  the  cell  can  produce  electric  current. 
However,  once  Crev  equals  Ece\\,  no  more  fuel  can  be  consumed 
and  no  more  current  can  be  produced.  Any  unused  fuels  in  the 
exhaust  reduce  the  fuel-cell  efficiency.  Fuel  utilization  £y  can 
be  written  as 


£y  =  1  — 


Wif.out  A/i  |  OL|t 
b/ 1'. in  A/if  ;n 


(6) 


where  the  “in”  and  “out”  refer  to  the  inlet  and  outlet  of  the  fuel 
cell.  The  Ah  refers  to  the  specific  enthalpy  associated  with  com¬ 
plete  oxidation  of  any  available  fuels.  This  definition  accounts 


Fig.  1.  Fuel-cell  efficiency  and  fuel  utilization  as  functions  of  operating  cell 
voltage  in  an  SOFC.  The  fuel  stream  is  based  on  the  equilibrium  state  of  50% 
H2  and  50%  CO  at  the  temperature  of  800  °C  and  atmospheric  pressure. 

for  the  energy  content  of  any  remaining  fuels  (or  fuel  byprod¬ 
ucts)  that  leave  in  the  fuel-cell  exhaust.  For  example,  even  though 
essentially  all  the  parent  fuel  (e.g.,  a  hydrocarbon)  is  consumed 
in  the  sense  it  is  no  longer  present  in  the  anode  exhaust,  there 
may  still  be  considerable  energy  available  in  the  form  of  other 
hydrocarbons  or  CO  and  H2 . 

Fig.  1  illustrates  the  typical  behavior  of  the  fuel-cell  efficiency 
and  fuel  utilization  as  functions  of  operating  cell  voltage  for  a 
given  fuel  stream  and  operating  conditions.  Operating  at  low 
cell  potential,  the  fuel  can  be  fully  utilized  and  the  fuel-cell 
efficiency  is  a  linear  function  of  cell  potential  £ceii.  The  slope  of 
this  function  depends  on  fuel  composition  as  well  as  operating 
temperature  and  pressure. 

The  reversible  efficiency  £r,  which  can  be  determined  ther¬ 
modynamically  (Eq.  (3)),  is  shown  as  a  circle  in  Fig.  1.  The 
reversible  efficiency,  which  would  be  achieved  by  operating  the 
cell  at  its  reversible  potential,  cannot  be  achieved  in  practice  be¬ 
cause  of  low  fuel  utilization  at  high  operating  potential.  The  ef¬ 
ficiency  reaches  a  maximum  at  an  operating  potential  of  around 
0.8  V.  As  discussed  below,  the  exact  conditions  for  maximum  ef¬ 
ficiency  depend  on  fuel  composition,  temperature,  and  pressure. 
The  quantitative  behavior  of  the  maximum  efficiency  versus  cell 
potential,  as  illustrated  in  Fig.  1 ,  was  explained  and  discussed 
by  Sidwell  and  Coors  using  a  somewhat  different  model  [7]  . 

2.  Thermodynamic  efficiency  model 

With  some  reasonable  assumptions,  the  efficiency  (i.e.,  Eq. 
(2))  may  be  evaluated  entirely  thermodynamically.  As  a  result, 
efficiency  can  be  computed  very  easily  over  wide  ranges  of  op¬ 
erating  conditions.  Model  development  is  facilitated  by  thinking 
of  fuel  and  air  flow  along  channels  separated  by  an  MEA  struc¬ 
ture.  As  the  fuel  flows  along  the  anode  channel,  it  is  continuously 
diluted  with  products  of  the  electrochemical  reactions  (primar¬ 
ily  CO2  and  H2O).  Consequently  the  local  reversible  potential 
/frev  decreases  along  the  length  of  the  channel.  Electrochemical 
charge  transfer,  and  hence  power  generation,  requires  that  Eiev 
exceed  the  operating  cell  potential  £ceii,  which  is  assumed  to 
be  spatially  uniform  along  the  channel  length.  If  the  channel 
is  sufficiently  long,  the  local  Erev  will  eventually  decrease  to 
operating  potential  Tfceii-  At  this  point  the  current  density  van¬ 
ishes  and  no  more  electric  power  can  be  produced.  The  effective 
length  for  power  generation  Leff  is  the  channel  length  at  which 
Erev  =  Cedi  -  It  is  at  this  point  where  the  fuel  utilization  £u  and 
fuel-cell  efficiency  s  are  maximum.  Although  Leff  varies  as  a 
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function  of  MEA  structure  and  operating  conditions,  it  is  inter¬ 
esting  to  note  that  both  eu  and  e  are  functions  only  of  the  inlet 
fuel  composition — they  are  independent  of  MEA  structure  and 
fuel  flow  rate.  However,  if  the  channel  length  is  less  than  Leff, 
the  actual  fuel  utilization  eu  and  efficiency  e  do  depend  on  MEA 
structure  and  flow  rates. 


2.1.  Reversible  cell  potential 

The  thermodynamic  analysis  begins  with  a  global 
electrochemical-oxidation  reaction  as 


K  K 

^^2v'lkXk  +  J2vo.kXk.  (7) 

4=1  4=1 


The  seemingly  complicated  nomenclature  is  needed  to  retain 

generality  in  the  fuel  and  oxidizer  mixtures  [8].  The  chemical 

symbol  for  the  Ath  species  (which  may  participate  as  a  fuel,  an 

oxidizer,  or  both)  is  written  as  Xk-  The  mole  fractions  within 

the  inlet  fuel  and  oxidizer  mixtures  are  given  as  n f,*  and  «0,4- 

Stoichiometric  coefficients  for  the  fuel  and  oxidizer  mixtures 

are  given  as  v'f  and  v'0,  and  stoichiometric  coefficients  of  the  Ath 

product  species  in  the  fuel  and  air  channels  are  given  as  v” ,  and 

v" , . 
o,k 

The  stoichiometric  coefficients  are  easily  determined  by  bal¬ 
ancing  the  reaction  for  particular  fuel  and  oxidizer  streams.  Each 
species  Xk  has  a  primary  identity  as  a  fuel,  an  oxidizer,  a  prod¬ 
uct,  or  an  inert,  and  this  identity  plays  an  essential  role  in  bal¬ 
ancing  the  reaction  to  determine  stoichiometric  coefficients  and 
in  assigning  charge  transfer  [8].  The  reactant  and  product  stoi¬ 
chiometric  coefficients  v  are  determined  to  balance  the  reaction. 
Depending  on  the  ion  transport  in  the  electrolyte,  products  may 
be  formed  in  the  fuel  or  the  oxidizer  channels.  For  example,  in 
a  polymer  electrolyte  membrane  (PEM)  cell  the  product  H2O  is 
formed  in  the  oxidizer  channel,  while  for  an  SOFC  the  product 
H2O  is  formed  in  the  fuel  channel.  This  difference  is  accom¬ 
modated  by  specifying  the  use  of  either  v'fk  or  v"  k  to  designate 
where  the  product  is  formed.  The  global  reaction  can  be  written 
in  compact  form  as 


Based  on  the  chemical-potential  balance,  the  reversible  cell 
potential  can  be  calculated  from  the  Nernst  equation  as 


AG 


neF 


T  ^  /  \ 

=  E°  H - -  V'(Vf«f,4  +  v'0n0,k  ~  vk)  In  (  —  )  ,  (10) 

»cf  \poJ 


where  E°  is  the  ideal  Nernst  potential  at  standard  con¬ 
ditions  (po  —  1  atm),  pk  the  partial  pressure  of  the  Ath 
species,  R=  8.314  J  mol-1  K-1  the  universal  gas  constant, 
and  F=96,485.309Cmol_1  is  the  Faraday  constant.  The  ideal 
Nernst  potential  at  the  standard  conditions  is  given  as, 


E° 


A  G° 
ntF  ’ 


(11) 


where  AG°  is  the  change  in  standard- state  Gibbs  free  energy 
between  products  and  reactants  of  the  global  electrochemical 
reaction  (Eq.  (8)).  Specifically, 


A  G° 


K 

+  vono,k  -  v'k)p,°k, 

4=1 


(12) 


where  pu°k  is  the  standard-state  chemical  potential  of  the  Ath 
species.  The  standard-state  thermodynamic  properties  of  ideal 
gases  depend  only  on  temperature  T,  and  are  readily  available 
in  databases  such  as  in  Chemkin  [9].  Assuming  H2  is  the  only 
electrochemically  active  fuel  species,  the  reversible  cell  poten¬ 
tial  can  be  written  as, 


ArevJU 

with 

AGn,  = 


AG 


IF 


H2  RT  \n PI  h 0,21 


IF  1/2  ’ 

PH2,aPo2,c 


H2  —  F-H20,a  —  Mh2,, 


(13) 


(14) 


It  is  well  known  that  CO  is  also  electrochemically  active.  As¬ 
suming  charge  transfer  via  H2  alone  is  justified  by  the  fact  that 
catalytic  reaction  of  CO  with  H2O  to  form  H2  and  CO2  is  fast. 
Thus,  the  water-gas-shift  process  to  produce  H2  competes  fa¬ 
vorably  with  the  relatively  slow  CO  electrochemical  oxidation. 


2.2.  Reaction  voltage 


K  K 

5>4X4^$>£x*,  (8) 

4=1  4=1 

where  v'k  =  v'fnftk  +  v'0n0,k  and  vk  =  Vfk  +  v”  k. 

The  number  of  electrons  transferred  by  the  global  electro¬ 
chemical  reaction  is  determined  from  the  half-cell  reactions  as 

K  K 

«e  =  ^2  vf”f,4Zf,4  =  J!  vo"o,4Zo,4-  (9) 

4=1  4=1 

where  zt,k  and  za,k  are  the  charges  associated  with  each  fuel  and 
oxidizer  species. 


Assuming  that  the  global  reaction  achieves  equilibrium 
within  the  cell,  it  follows  from  species  mass  balances  in  fuel 
and  oxidizer  channels  that 


hlf,4,in  bl  f.UouL  — 


ni  0,4.  in  ril  0,4,  out  — 


Wk(v'fnik  - 

Tl^r 

Wk(v'0n0,k  -  v^k)-1^-. 

Yl^r 


(15) 

(16) 


The  net  electric  current  produced  within  the  cell  is  /e .  The  species 
mass  flow  rates  at  the  inlets  and  outlets  of  the  fuel  and  oxidizer 
channels  are  represented  as  riik,  and  Wk  are  species  molecu¬ 
lar  weights.  Multiplying  both  of  these  equations  by  the  species 
specific  enthalpies  (i.e.,  enthalpy  per  unit  mass),  which  for  a 
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uniform- temperature  system  are  constants,  and  summing  over 
all  species  yields 

h  K 

0 mh)fM  -  (m/0f,out  =  E  Wkhk(v'fntk  -  v" k),  (17) 

ntt  k= l 

h  K 

(mh) 0,in  -  (m/x)0  0ut  =  — ^  ^2  Wkhk(v'0no,k  -  v”  k). 

k=  1 

(18) 


complete  oxidation  of  the  inlet  fuel  stream  can  be  represented 
in  terms  of  the  enthalpies  of  fuel,  oxidizer,  and  final  products  as 

(m  A/t)f, i„  =  (mh) f,in  +  (mA)*£  -  (m/t)^,  (22) 

where  (m h)^  is  an  “enthalpy  rate”  of  the  oxidizer  required  to 
fully  oxidize  the  inlet  fuel  stream  and  (mh)^  is  the  enthalpy  rate 
of  the  products  of  the  completely  oxidized  initial  fuel  stream. 
An  analogous  relationship  applies  to  the  exhaust  stream,  which 
may  contain  unspent  fuel, 

(m  Ah)um  =  (mh) f.out  +  (mh)™nt  -  (mh)™nV  (23) 


The  total  mass  flow  rates  and  enthalpies  in  each  channel  are 
represented  as  m  =  Y2  rhk  and  h  —  Y2  Ykhk,  where  Yk  =  mk/m 
are  the  species  mass  fractions.  Adding  these  two  equations  and 
multiplying  by  the  cell  operating  potential  /;cc]i  yields 


Ece\\[(mh)m  -  (mh) out]  =  E^We, 


(19) 


where  (rhh)m  —  (mh) f,;n  +  (mh\ ;n  and  (mh)out  —  (mh) f,out  + 
(mh) 0,out-  The  net  electric  work  (power)  is  written  as  We  — 
/cTccii  and  a  “reaction  voltage”  is  defined  as 


k=  1 


Wkhk 

rie.F 


[(v’fif'k  +  v'0n0,k)  -  (v'!k  +  u"  j.)] 


1 

neF 


K 


E  Wkhk[v'k  -  Vk]  =  ~ 
k=  1 


AH 

neF 


(20) 


2.3.  Thermodynamic  efficiency 

With  the  reversible  potential  written  as  Elev  =  —  A G/neF 
and  using  Eq.  (19),  the  net  cell  efficiency  can  be  written  as 

e  _  We  _  f  £rey)  (  Ecell  \  (  (mh)in  -  (mh)oul  \ 
ihf.in  A/z f,in  \  /-),  /  \  /-rev  /  \  ffl f  in  All f  jn  J 

=  SRfiveu-  (21) 

In  other  words,  the  cell  efficiency  is  the  product  of  a  reversible 
efficiency  Sr,  a  voltage  efficiency  f:\.  and  a  utilization  efficiency 

f?U- 

Eq.  (21)  is  a  closed-form  expression  for  fuel-cell  efficiency. 
However,  as  discussed  below,  using  Eq.  (21)  requires  evaluat¬ 
ing  (mh)out,  which  requires  assuming  that  the  exhaust  stream 
is  in  chemical  equilibrium.  It  is  important  to  note  that  Eq.  (21) 
does  not  depend  on  cell  structure  or  internal  polarization  losses. 
It  does  account  for  the  fuel  mixture,  operating  temperature  and 
pressure,  fuel  utilization,  and  cell  voltage.  Achieving  the  ther¬ 
modynamically  predicted  efficiency  requires  that  the  cell  be  as 
large  as  needed,  with  larger  areas  needed  for  cells  with  higher 
internal  polarization  losses. 

2.4.  Equivalent  expressions  for  utilization 

The  definition  for  eu  in  Eq.  (21)  appears  to  be  different  from 
the  one  in  Eq.  (6).  However,  some  further  manipulation  reveals 
that  they  are  indeed  the  same.  The  potential  heat  release  from 


For  a  given  fuel  stream,  the  fully  oxidized  products  must  have 
the  same  enthalpy  rates  at  the  cell  inlet  and  the  cell  out¬ 
let  (i.e.,  (m/z)pUjn  =  (m/t)pUQUt).  Furthermore,  the  difference  in 
the  oxidizer  required  to  completely  oxidize  the  inlet  and  out¬ 
let  streams  must  be  provided  from  the  cathode  stream  (i.e., 
(mh)^  -  (»i/z)£oUt  =  (rhhfo.in  -  (mh) 0,out)-  It  follows  that 

(m  A/i)0,in  -  (m  A/z)f,out 

=  (mh) 0,in  -  (mh) f,out  +  (rhh)0;m  -  (mh) 0,0ut 

=  (mh)in  -  (mh)oat.  (24) 

Thus,  the  two  representations  of  eu  in  Eqs.  (6)  and  (21)  are  the 
same. 


2.5.  Algorithm  for  efficiency  calculation 


With  two  assumptions,  the  maximum  fuel  utilization  can  be 
determined.  Assume  first  that  the  fuel  cell  is  sufficiently  large 
that  all  the  fuel  that  can  be  consumed  is  consumed.  In  other 
words,  the  fuel  content  of  the  exhaust  stream  has  been  depleted 
and  diluted  to  the  point  that  ETev  =  ECeii-  Assume  further  that 
the  exhaust  mixture  is  in  chemical  equilibrium.  For  a  given  £ceii> 
the  exhaust  composition  can  be  determined  via  the  relationship 
between  composition  and  Elev.  Because  it  is  difficult  to  evalu¬ 
ate  directly  the  numerator  in  Eq.  (6),  the  following  approach  is 
followed. 

The  species  molar  flow  rates  at  the  inlet  to  the  fuel  and  ox¬ 
idizer  channels  are  specified  as  Nfnk  and  N1fk.  Assuming  that 
some  fraction  y  of  the  inlet  fuel  is  electrochemically  oxidized 
via  the  global  reaction  (Eq.  (7)),  the  outlet  molar  flow  rates  may 
be  determined  as 


liroul  _ 

,k  — 


N ] 


f,  tot 


[(1  -  y)v'fn{'k  +  yvfk]. 


Nln 

ArOUt  A/in  |  f,  tOt  r  Jt  ,  /  „  n 

No.k  —  No.k  H  77  Y(-Vo.k  -  WtrL 


(25) 

(26) 


where  /Vp"ot  is  the  total  molar  flow  of  the  fuel  stream  at  the  inlet. 
In  this  way,  both  the  fuel  and  air  flow  rates  at  the  outlet  can  be 
defined  in  terms  of  a  single  parameter  y. 

To  calculate  the  reversible  potential  Erev,  both  the  outlet  fuel 
stream  and  air  stream  are  assumed  to  be  in  chemical  equilibrium. 
The  equilibrium  species  distributions  are  computed  assuming 
the  element  balances  associated  with  the  molar  flow  rates  at  the 
channel  outlets,  Nff  and  N°uk.  Thus,  once  the  inlet  flows  and  the 
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global  reaction  is  specified,  Elev  can  be  determined  in  terms  of 
y  alone.  By  assuming  that  charge  transfer  proceeds  entirely  via 
H2  oxidation,  E[ev  is  determined  from  Eq.  (13)  with  the  partial 
pressures  taken  from  the  equilibrium  compositions  of  the  anode 
and  cathode  outlets. 

Once  the  reversible  potential  Elev  equals  the  cell  operating 
potential  Eceu,  no  further  fuel  can  be  consumed  electrochem- 
ically  and  no  further  power  is  generated.  The  maximum  fuel 
consumption  ymax  can  be  determined  from  the  iterative  solution 
of 

Erev(Kmax)  —  Ece!l  •  (27) 

Once  ymax  is  determined,  the  corresponding  compositions  of  the 
fuel  and  air  streams  at  the  outlets  can  be  calculated.  With  the  fuel 
composition  at  the  outlet  known,  the  heating  value  associated 
with  full  oxidation  of  the  remaining  fuel  can  be  computed.  Once 
ih f)0ut  A/if  out  is  computed,  ey  is  determined  from  Eq.  (6).  The 
net  efficiency  then  follows  from  Eq.  (21). 


2.6.  Effect  of  mass  flow  rate  on  efficiency 


Assuming  that  the  cathode-side  air  flow  rate  is  sufficiently 
high  that  maximum  fuel  utilization  can  be  achieved,  the  fuel¬ 
cell  efficiency  can  be  shown  to  be  independent  of  the  flow  rate 
of  the  inlet  fuel  stream.  The  mass  flow  rate  at  the  fuel-channel 
outlet  can  be  represented  as, 

K 

flif.out  =  NffWk,  (28) 

k=  l 

Based  on  Eq.  (25),  the  mass  flow  rate  can  be  rewritten  as, 

K  w 

w'if.out  =  Nftot  -T I K] 1  -  y>'fn  u  +  Wfk\'  (29) 

k=  1  vf 
or 


where  Wf;ia  =  J2k=i  nf,kW/c  is  the  mean  molecular  weight  of 
the  inlet  fuel  mixture.  Therefore,  the  ratio  of  the  mass  flow  rates 
at  the  outlet  and  the  inlet  of  the  fuel  channel  can  be  written  as 

™f,out  ,,  ,  ,  Ef=l(Vf'jt/Vf)WJfc 

- —  =  (i  -y)  +  y  — K  - ■  (31) 

mi,in  J2k=  1  nf,kWk 

It  is  apparent  the  ratio  of  the  mass  flow  rates,  Wf,outA«f,in,  de¬ 
pends  only  on  the  composition  of  the  inlet  fuel  mixture,  the 
global  electrochemical  reaction,  and  the  fuel  utilization;  it  does 
not  depend  on  the  flow  rate  through  the  fuel  channel.  Thus,  at 
/max  where  Eiev  =  Ece n,  Eqs.  (6),  (21),  and  (31)  indicate  that 
both  the  fuel-cell  efficiency  and  the  fuel  utilization  are  indepen¬ 
dent  of  the  mass  flow  rates  through  the  fuel  and  air  channels. 


2.7.  Effects  of  cell  voltage  and  fuel  stream  on  efficiency 

Fig.  2  shows  maximum  efficiencies  and  fuel  utilization  as 
functions  of  the  operating  cell  voltage  for  several  fuel  streams. 


Fig.  2.  The  predicted  maximum  thermodynamic  fuel-cell  efficiency  and  fuel 
utilization  as  functions  of  cell  voltage  for  different  fuel  streams  at  800  °C  and 
1  atm.  The  fuel  streams  illustrated  here  include  the  pure  H2  fuel  stream,  the 
equilibrium  composite  fuel  stream  based  on  CH4-steam  and  JP8-steam  streams 
with  the  steam-carbon  ratio  of  1  (CH4-I,  JP8-1)  and  3  (CH4-3,  JP8-3)  at  800  °C 
and  1  atm.  The  lower  panel  is  an  expanded  view  for  the  operating  voltages  from 
0.7  V  to  IV. 


These  include  pure  H2  and  mixtures  of  CH4-steam,  and  JP8- 
steam.  Steam-carbon  ratios  of  1  and  3  are  considered,  designated 
as  CH4-I,  CH4-3,  JP8-1,  and  JP8-3.  JP8  is  assumed  to  be  a 
mixture  of  80.1%  Ci3H28,  1.9%  Ci3H26,  and  18%  Ci3Hi4.  For 
the  hydrocarbon-steam  mixtures,  the  fuel  that  enters  the  fuel 
cell  is  the  equilibrium  composition  at  the  fuel-cell  temperature 
and  pressure.  Table  1  shows  the  equilibrium  major-species  mole 
fractions  at  800  °C  and  atmospheric  pressure  for  the  four  fuel- 
steam  mixtures. 

In  all  cases  shown  in  Fig.  2,  the  cell  is  assumed  to  be 
isothermal  at  800  °C,  atmospheric  pressure,  and  the  cathode 
stream  is  undiluted  air.  The  effects  of  the  three  contributing 
efficiencies  are  evident  in  the  figure.  At  low  cell  voltage  the 
efficiency  increases  linearly  as  a  result  of  e\ .  In  this  region 
the  utilization  is  nearly  100%  and  the  slope  of  the  efficiency 
function  depends  on  the  reversible  efficiency  £r.  At  around 
0.8  V  the  utilization  begins  to  fall,  contributing  to  a  decreas¬ 
ing  efficiency.  The  utilization  decreases  because  the  operat¬ 
ing  voltage  exceeds  the  reversible  potential  of  the  partially 
depleted  fuel  stream.  At  sufficiently  high  voltage,  utilization 
(and  hence  efficiency)  decreases  to  zero  when  none  of  the 
fuel  can  be  electrochemically  oxidized.  The  cell  efficiency  and 


Table  1 

Gas-phase  equilibrium  mole  fractions  for  four  fuel-steam  mixtures  at  800  °C 
and  1  atm 


CH4-1 

CH4-3 

JP8-1 

JP8-3 

ch4 

0.02674 

0.00042 

0.02764 

0.00026 

CO 

0.22968 

0.10792 

0.31686 

0.11644 

H2 

0.71684 

0.55761 

0.62793 

0.48120 

h2o 

0.01979 

0.27558 

0.01775 

0.31781 

co2 

0.00695 

0.05846 

0.00982 

0.08430 
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Fig.  3.  Illustration  of  how  the  fuel-cell  efficiency  varies  as  a  function  of  the 
operating  cell  potential  for  a  given  fuel  utilization. 

fuel  utilization  of  CH4-steam  and  JP8-steam  fuel  streams  with 
the  steam-carbon  ratio  of  1  are  almost  the  same  as  those  for 
the  pure  H2  fuel  stream.  However,  as  more  steam  is  added 
to  the  hydrocarbon  fuels,  the  cell  efficiency  drops  slightly, 
and  the  operating  voltage  for  maximum  efficiency  also  drops 
slightly. 

There  are  often  good  reasons  to  operate  the  fuel  cell  itself  at 
less  than  maximum  efficiency,  seeking  to  optimize  overall  effi¬ 
ciency  of  a  full  system.  For  example,  excess  fuel  in  the  fuel-cell 
exhaust  may  be  burned  to  provide  heat  needed  in  a  recuperator 
to  preheat  the  air  stream  prior  to  entering  the  fuel  cell.  Thus,  the 
maximum  efficiency  and  fuel  utilization  shown  in  Fig.  2  may  not 
represent  the  best  operating  point.  Rather,  operating  at  reduced 
fuel  utilization  may  be  preferred.  Fig.  3  illustrates  the  relation¬ 
ships  of  the  maximum  efficiency  and  utilization  to  efficiency  at 
a  specified  utilization. 

For  fixed  utilization  (illustrated  in  Fig.  3  as  70%),  the 
maximum-utilization  function  can  be  used  to  determine  the  max¬ 
imum  operating  potential  above  which  the  fixed  fuel  utilization 
can  no  longer  be  achieved.  At  this  maximum  operating  potential, 
the  maximum-efficiency  curve  shows  the  maximum  efficiency 
for  the  specified  utilization.  The  straight  line  connecting  the  ori¬ 
gin  to  this  maximum  efficiency  represents  the  efficiency  at  fixed 
utilization. 

At  a  given  utilization,  it  is  useful  to  know  the  maximum  ef¬ 
ficiency  and  operating  potential  needed  to  achieve  the  maxi¬ 
mum  efficiency.  Fig.  4  shows  maximum  efficiency  as  a  function 
of  utilization  for  two  methane-steam  fuel  mixtures.  The  figure 
also  shows  the  operating  potential  needed  to  achieve  the  max¬ 
imum  efficiency.  The  information  in  Fig.  4  is  identical  to  that 
in  Fig.  2;  it  simply  re-plots  the  data  using  utilization  as  the  ab¬ 
scissa. 


Fig.  4.  Maximum  fuel-cell  efficiency  and  the  corresponding  operating  poten¬ 
tial  as  functions  of  maximum  fuel  utilization  for  CFLi-steam  equilibrium  fuel 
streams  at  a  steam-carbon  ratios  of  1  (CH4-I)  and  3  (CH4-3).  In  all  cases  the 
cell  is  at  1  atm  and  800  °C. 


CH4  /  (CH4  +  H2O) 


Fig.  5.  Efficiency  characteristics  for  SOFCs  with  methane-steam  fuel  mixtures. 
The  abscissa  is  the  mole  fraction  of  methane  in  methane-steam  mixture.  The 
fuel  that  enters  the  fuel  cell  itself  is  the  equilibrium  composition  that  results 
from  the  given  methane-steam  mixture  at  the  fuel-cell  temperature  and  atmo¬ 
spheric  pressure.  The  top  panel  shows  maximum  possible  efficiency.  The  middle 
panel  shows  the  utilization  at  maximum  efficiency.  The  lower  panel  shows  the 
operating  potential  needed  to  achieve  maximum  efficiency.  The  solid  markers 
represent  the  points  at  which  equilibrium  predicts  solid-carbon  formation. 

3.  Efficiency  maps 

Because  the  thermodynamic  models  are  so  efficient  compu¬ 
tationally,  it  is  easy  to  map  large  regions  of  parameter  space. 
This  can  be  helpful  in  identifying  trends,  seeking  optimal  oper¬ 
ating  conditions.  Fig.  5  shows  SOFC  maximum-efficiency  char¬ 
acteristics  as  functions  of  methane-steam  mixtures  and  cell  op¬ 
erating  temperature.  Clearly  maximum  efficiency  increases  as 
methane  content  increases.  Reducing  temperature  also  increases 
efficiency  significantly.  In  all  cases  utilization  is  above  90%  at 
maximum  efficiency.  The  operating  potential  to  achieve  maxi¬ 
mum  efficiency  increases  substantially  as  operating  temperature 
decreases.  The  solid  markers  on  the  curves  represent  the  point 
at  which  chemical  equilibrium  predicts  the  formation  of  solid 
carbon.  At  methane  concentrations  above  these  levels,  the  sys¬ 
tems  design  needs  to  consider  coke  formation  that  may  degrade 
cell  performance.  If  higher  hydrocarbons  also  enter  the  fuel  cell 
predicting  deposit  formation  is  more  complex.  In  addition  to 
surface-catalyzed  coke  formation,  gas-phase  routes  can  lead  to 
polyaromatic  hydrocarbon  deposits  [10,11], 

4.  Channel  flow  model  and  power  density 

It  is  interesting  to  compare  the  results  of  the  thermodynamic- 
efficiency  model  with  a  much  more  detailed  model  that  directly 
incorporates  electrochemical  charge-transfer  kinetics,  heteroge¬ 
neous  reforming  kinetics,  fuel  and  air  flow  in  anode  and  cath- 
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ode  channels,  dense-electrolyte  ion  transport,  and  porous  media 
transport  in  electrodes  [4].  The  detailed  model  confirms  the  pre¬ 
dictions  of  the  thermodynamic  model.  The  detailed  model  also 
provides  much  more  information,  such  as  power  density  and 
thermal  loads.  The  detailed  model  also  provides  a  great  deal  of 
local  information  throughout  the  cell,  none  of  which  the  thermo¬ 
dynamic  model  can  deliver.  Unlike  the  thermodynamic  model, 
the  detailed  model  requires  very  specific  information  about  the 
cell  architecture  and  the  ME  A,  as  well  as  operating  conditions. 
For  the  sake  of  illustration  the  cell  described  by  Zhu  et  al.  [4]  is 
used  here  in  the  detailed  model.  Summarized  briefly,  the  anode  is 
a  1.22-cm-thick  porous  Ni-YSZ  cermet.  The  dense  electrolyte  is 
8-YSZ  that  is  25  |xm  thick.  The  cathode  is  a  30  |xm  porous  LSM 
structure.  For  the  example  here,  the  detailed  model  uses  square 
channel  cross-sections  of  1  mm2.  The  electrochemically  active 
perimeter  is  1.2  mm,  which  is  0.2  mm  wider  than  the  channel 
electrode  interface.  Other  physical  and  chemical  parameters,  as 
well  as  details  of  the  model,  can  be  found  in  [4], 

The  upper  panel  in  Fig.  6  illustrates  the  performance  of  the 
nominal  MEA  structure  operating  in  a  humidified  H2  fuel  stream 
and  a  diluted  H2  fuel  stream  at  800  °C  and  1  atm.  To  assist  in¬ 
vestigating  the  effects  of  the  MEA  performance  on  the  fuel-cell 
efficiency,  a  reduced-performance  MEA  was  constructed.  Leav¬ 
ing  all  other  cell  parameters  unchanged,  the  exchange  current 
densities  at  the  anode-electrode  interface  ijj^  and  at  the  cathode- 
electrolyte  interface  i^n  are  both  reduced  to  1  A  cm-2  (see  Zhu 
et  al.  [4]  for  detailed  definitions).  The  lower  panel  in  Fig.  6  shows 
the  MEA  performance  is  significantly  degraded. 

Fig.  7  compares  directly  results  from  the  thermodynamic  and 
detailed  models  for  a  cell  operating  on  a  fuel  stream  mixture 
of  12%  CH4,  66%  H2,  and  22%  CO,  which  is  essentially  the 
gas-phase  chemical  equilibrium  composition  of  an  initial  mix- 


Fig.  6.  Polarization  characteristics  for  two  MEA  structures,  each  operating  on 
H2-H2O  fuel  mixtures.  The  upper  panel  shows  performance  of  the  MEA  that  is 
described  in  Zhu  et  al.  [4].  The  lower  panel  shows  a  lower-performance  MEA. 
In  both  cases  the  cells  are  operating  at  800  °C  and  atmospheric  pressure. 


Fig.  7.  Comparison  of  the  maximum  fuel-cell  efficiency  and  fuel  utilization 
calculated  from  the  thermodynamic  model  and  a  detailed  model  [4].  The  inlet 
fuel  stream  is  the  equilibrium  composition  resulting  from  an  initial  fuel  mixture 
of  60%  CH4  and  40%  H2O  at  800  °C  and  1  atm.  The  fuel  cell  is  operated  at 
800  °C  and  1  atm.  The  overall  power  densities  for  the  normal  MEA  and  the 
performance-reduced  MEA  are  calculated  with  the  detailed  model. 

ture  of  60%  CH4  and  40%  FFO  at  800  °C  and  1  atm.  The  inlet 
velocity  is  10  cm  s~  1 .  The  models  predict  virtually  identical  ef¬ 
ficiency  and  utilization  performance.  This  excellent  comparison 
is  observed  over  very  large  ranges  of  fuels  and  operating  con¬ 
ditions.  Furthermore,  the  efficiency  and  utilization  predicted  by 
the  detailed  model  are  virtually  identical  with  either  the  nom¬ 
inal  or  reduced-performance  MEA.  However,  as  discussed  be¬ 
low,  when  the  reduced-performance  MEA  is  used  the  cell  must 
be  much  larger  to  achieve  the  maximum  efficiency.  Fig.  7  also 
compares  the  power  density  at  the  maximum  fuel  utilization  for 
the  two  MEA  structures.  As  expected,  the  power  density  for 
the  performance-reduced  MEA  is  much  lower  than  that  for  the 
normal  MEA.  The  different  power  densities  lead  to  different 
channel  lengths  to  achieve  the  same  maximum  efficiency.  For 
the  conditions  in  these  examples,  the  nominal  MEA  results  in  a 
channel  length  of  9. 1  cm  to  achieve  the  maximum  efficiency  of 
63%.  The  reduced-performance  MEA  requires  a  much  longer 
channel  length  of  17.9  cm. 

It  is  interesting  to  note  that  the  power  densities  show  max¬ 
ima  as  functions  of  operating  cell  potential.  The  voltage  for 
the  peak  power  density  is  around  0.68  V  for  the  nominal  MEA 
in  this  example,  and  is  reduced  slightly  to  about  0.63  V  for 
the  performance-reduced  MEA.  In  both  cases  the  peak  power 
density  occurs  at  operating  potentials  well  below  the  operat¬ 
ing  potential  for  maximum  efficiency.  From  the  fuel-cell  de¬ 
sign  and  optimization  point  of  view,  the  fuel  cell  should  be  de¬ 
signed  and  operated  to  achieve  appropriate  tradeoffs  between 
efficiency,  utilization,  and  power  density.  The  models  described 
here  provide  tools  that  can  be  very  effective  in  evaluating  these 
tradeoffs. 

5.  Conclusions 

Two  models  are  developed  and  used  to  predict  SOFC  effi¬ 
ciencies  for  a  wide  range  of  fuels  and  operating  conditions.  The 
net  efficiency  is  the  product  of  a  reversible  efficiency,  voltage 
efficiency,  and  fuel  utilization.  Cell  efficiency  is  a  strong  func¬ 
tion  of  operating  voltage,  with  maximum  efficiency  achieved  at 
cell  voltages  in  the  range  of  0.8  V.  The  thermodynamic  model 
predicts  best  possible  performance,  independent  of  specific  cell 
design.  The  channel-based  model  incorporates  details  of  MEA 
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structure,  fluid  transport,  and  charge-transfer,  and  reforming  ki¬ 
netics.  Both  models  can  be  used  to  guide  the  optimization  of  cell 
design  and  operation. 
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